The First High-Precision Radial Velocity Search for 

Extra-Solar Planets 



Gordon A. H. Walker 

1234, Hewlett Place, Victoria BC, V8S 4P7, Canada 



Abstract 

The reflex motion of a star induced by a planetary companion is too small to 
detect by photographic astrometry. The apparent discovery in the 1960s of 
planetary systems around certain nearby stars, in particular Barnard's star, 
turned out to be spurious. Conventional stellar radial velocities determined 
from photographic spectra at that time were also too inaccurate to detect 
the expected reflex velocity changes. In the late 1970s and early 1980s, the 
introduction of solid-state, signal-generating detectors and absorption cells to 
impose wavelength fiducials directly on the starlight, reduced radial velocity 
errors to the point where such a search became feasible. Beginning in 1980, 
our team from UBC introduced an absorption cell of hydrogen fluoride gas 
in front of the CFHT coude spectrograph and, for 12 years, monitored the 
radial velocities of some 29 solar-type stars. Since it was assumed that extra- 
solar planets would most likely resemble Jupiter in mass and orbit, we were 
awarded only three or four two-night observing runs each year. Our survey 
highlighted three potential planet hosting stars, 7 Cep (Kl IV), f3 Gem (K0 
III), and e Eri (K2 V). The putative planets all resembled Jovian systems 
with periods and masses of: 2.5 yr and 1.4 Mj, 1.6 yr and 2.6 Mj, and 6.9 
yr and 0.9 Mj, respectively. All three were subsequently confirmed from 
more extensive data by the Texas group led by Cochran and Hatzes who also 
derived the currently accepted orbital elements. 

None of these three systems is simple. All 5 giant stars and the supergiant 
in our survey proved to be intrinsic velocity variables. When we first drew 
attention to a possible planetary companion to 7 Cep in 1988 it was classified 
as a giant, and there was the possibility that its radial velocity variations and 



Email addresses: gordonwaOuvic . ca (Gordon A. H. Walker) 



Preprint submitted to Elsevier 



December 16, 2008 



those of (3 Gem (KO III) were intrinsic to the stars. A further complication 
for 7 Cep was the presence of an unseen secondary star in an orbit with 
a period initially estimated at some 30 yr. The implication was that the 
planetary orbit might not be stable, and a Jovian planet surviving so close 
to a giant then seemed improbable. Later observations by others showed the 
stellar binary period was closer to 67 yr, the primary was only a sub-giant 
and a weak, apparently synchronous chromospheric variation disappeared. 
Chromospheric activity was considered important because k 1 Cet, one of our 
program stars, showed a significant correlation of its radial velocity curve 
with chromospheric activity. 

e Eri is a young, magnetically active star with spots making it a noisy 
target for radial velocities. While the signature of a highly elliptical orbit 
(e = 0.6) has persisted for more than 3 planetary orbits, some feel that even 
more extensive coverage is needed to confirm the identification despite an 
apparent complementary astrometric acceleration detected with the Hubble 
Space Telescope. 

We confined our initial analyses of the program stars to looking for cir- 
cular orbits. In retrospect, it appears that some 10% of our sample did in 
fact have Jovian planetary companions in orbits with periods of years. 

Key words: stars: activity, stars: individual: 7 Cep, (3 Gem, e Eri,, stars: 
late- type, stars: starspots, stars: rotation, stars: chromospheres,, stars: 
exoplanets 



1. The very idea of searching for extra-solar planets in the 1970s 

It is quite hard nowadays to realise the atmosphere of skepticism and 
indifference in the 1980s to proposed searches for extra-solar planets. Some 
people felt that such an undertaking was not even a legitimate part of as- 
tronomy. It was against such a background that we began our precise radial 
velocity (PRV) survey of certain bright solar-type stars in 1980 at the Canada 
France Hawaii 3.6-m Telescope (CFHT). Up to that time, the only serious 
searches for extra-solar planets had been astrometric studies that looked for 
long term periodic displacements of stars on the sky caused by unseen com- 
panions. Unfortunately, the tiny angular displacements were much smaller 
than the proper motions and parallaxes for nearby stars and beyond the 
capability of ground-based astrometric telescopes using photographic emul- 
sions. 
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For a circular planetary orbit, the semi-angular reflex displacement, 9, of 
the parent star on the sky caused by a companion planet is: 



Q ~ — EL milli-arcsec, (1) 

D M st v ' 

where R p i is the distance of the planet from its parent star in AU, M v \ and 
M st are the planet and stellar masses in Jovian and solar masses, respectively. 
D is the distance of the star from the Earth in parsecs. 6 increases linearly 
with R p i (and as the | power of the period) . For the Sun- Jupiter system seen 
from a distance of 5 pc, 9 would be ~1 milli-arcsec and take nearly 12 years to 
complete a cycle. Considering that the average ground-based photographic 
image was spread over >1000 milli-arcsec, the challenge of detecting such a 
miniscule displacement over several decades is obvious. 

Nonetheless, tentative extra-solar planetary discoveries were claimed for 
several nearby stars from astrometry with the Sproul Observatory 24 inch 
refractor, but the signals were close to the error s of measurement. T he best 



known identifications were for Barnard's Star (Ivan de Kampl . I1969T). The 



orbita l periods in all cases were close to multiples of 8 years, and iHershey 



( 119731 ) demonstrated that the epochs of the 'discontinuities' in the Sproul 
data coincided with times when the telescope primary lens had been removed 
for cleaning. 

The amplitude of the periodic reflex variation of a planet-hosting star's 
radial velocity complements the astrometric variation. Again, for circular 
orbits, the amplitude of the stellar radial velocity change, 5RV, would be: 

30M pl sim x 

dHV ~ (M st R pl y/* km ' {2) 

where % is the inclination of the orbital axis to the line of sight and the 
units are the same as those in equation 1. The radial velocity excursion is 
independent of the distance from the Earth and increases inversely as \JR p i- 
As a star, 5RV for the Sun as caused by Jupiter would be ~ ±13 m s _1 over 
12 years when viewed at 2=0. Unlike astrometry, radial velocities favor the 
detection of tightly bound planets. 

In the early 1970s, quality radial velocities determined from photographic 
plates had typical errors of ~1 km s _1 . The large error was in part due to 
the photographic process but, more seriously, to the effect of guiding errors 
at the spectrograph slit which cause the stellar spectrum to shift relative to 
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the bracketting comparison arc lines. Obviously, an improvement of more 
than an order of magnitude coupled with a long term stability of years, was 
necessary if radial velocities were to be useful in looking for extra-solar Jovian 
analogs. 

The above comments also reflect the belief at the time that planetary 
systems beyond the solar system would most likely be similar to our own, a 
belief strongly bolstered by theoretical simulations even down to the Titus 
Bode relation (jlsaacman fc Saganl . 119771 )! No one predicted the populations 
of Jovian planets in few-day, or highly elliptical, orbits which eventually 
turned up. As a result, early search strategies concentrated on long-term 
monitoring with observations spaced out over months and years. 



2. The UBC PRV Planet Search 

Development of digital TV systems for astronomical spectroscopy began 
at UBC in 1969. These evolved from conventional cameras with electron- 
beam read-out to self-scanned silicon diode arrays (I Walkerl . 119771 ). Access 
to the horizontal McKellar coude spectrograph fed by the 1.2-m Dominion 



Astrophysical Observatory (DAO) telescope 



which to test our cameras ([Richardson et al. 



j rovid ed a perfect laboratory in 



1972). 



With conventional TV cameras, electron reading-beam jitter and beam- 
bending are serious problems in maintaining fidelity of stellar line positions 
in the digitised spectra. These errors were in addition to those from guiding 
- although the introduction of the Richardson (1966) image sheer s ignifi- 



cantly reduced the latter. Glaspey and Fahlman (see I Walker et all Il973l ) 
demonstrated that one could use telluric water vapor lines to correct for all 
of these shifts, raising the possibility of achieving the true radial velocity 
precision inhe rent in the origina l data . The use of telluric lines was not a 



new idea (e.g. I Griffin &: Grifhnl . Il973l ). but it was particularly appropriate 



to the case of digital spectra of wide dynamic range and high signal to noise, 
and I realised that it raised the possibility of sufficient precision to look for 
extra-solar planets. 

The introduction of self-scanned diode arrays eliminated all of the prob- 
lems with electron beams and made the planet search a serious possibility. 
Bruce Campbell who joined me as a Post-Doctoral Fellow picked up on the 
idea and took it a step further by suggesting a captive gas. On the rec- 
ommendation of Alexander Douglas and Gerhard Herzberg of the Canadian 
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National Research Council Herzberg Institute of Astrophysics, hydrogen flu- 
oride (HF) was chosen because the 3-0 vibration band has a well spaced 
comb of lines with no confusion from isotopic shifts, the lines are of similar 
natural width to those in the stellar spectra, and there are no telluric lines 
contaminating th e A8700 region. 

Bruce and I (jCampbell &: Walkerl . 1 19791 ) set up an absorption cell just 
ahead of the McKellar spectrograph slit on 27 December 1978 to try the first 
'stellar' observations. Frankly, it was quite unsafe. HF is highly corrosive 
and toxic. The cell had to be heated to 100 C to prevent the HF from 
polymerising and the cell windows being plexiglass (HF attacks glass) warped 
with the heat. Nonetheless, we took a series of exposures of the Sun with the 
telescope mirror covers closed - enough light got through the gaps between 
the covers to give us a good signal from the 1024 Reticon diode array. 

It worked! Thanks to the small spectral shift program developed by 
Fahlman and Glaspey we were able quickly to analyse the spectral series and 
demonstrate that HF fiducial lines could in principle reduce radial velocity 
errors by the necessary order of magnitude from those being determined pho- 
tographically. One important aspect of the captive gas technique (as opposed 
to the use of telluric lines) is that stellar spectra can also be taken with the 
absorption cell out of the light beam to provide a template spectrum against 
which to determine differential velocities. Further, unlike the situation with 
telluric lines, spectra of the captive gas can be taken independently with a 
continuum source to allow the HF lines to be carefully aligned and removed 
from the stellar spectra. Exam ples of stellar spectra with and without the 



HF lines ar e show n in Figure 1 (jCampbell et all 119881 ) 



Beckers! (119761 ) was probably the first to use a captive g as, I2, to deter- 



mine wavelength shifts accur ately in solar spectra (see also iKoch &: Wohll . 



19841 ) . iMarcy fc Butlerl (119921 ) later adopted I2 for their very extensive PRV 



program. 

It should be emphasised that our radial velocities were precise but not 
accurate, and only appropriate for measuring differential velocities or the 
accelerations of the target stars. The wavelengths of stellar absorption lines 
depend on conditions in the stellar atmosphere in addition to the star's actual 
motion. There are upwelling velocities ~1 km s _1 of gasses hotter than 
their surroundings as well as cooler descending material, leading to the well 
known granulation pattern on the Sun. We had to assume that the average, 
weighted surface velocity remained constant (a doubtful assumption in many 
cases). The Earth itself is in revolution at 30 km s _1 about the solar system 
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barycentre and measured velocities must be corrected for this effect, which 
can amount to tens of km s _1 between observing runs and changes of several 
m s _1 during a single exposure. We had to accurately correct for the latter 
by estimating the mean time of the exposure weighted by the flux rate of 
starlight reaching the detector. 

The wavelength region covered by the HF lines (Figure 1) included one 
of the red Ca II triplet lines at A8662. The strength of this line is sensi- 
tive to chromospheric activity. Such activity introduces reversals of variable 
intensity within the Ca II lines. We considered this an important bonus 
because stellar granulation patterns might be affected by magnetic cycles 
and starspots thereby introducing apparent changes in radial velocity and 
we would be in a position to calibrate the effect if it existed. The period of 
Jupiter's orbit and the sunspot cycle are similar for instance. 

Bruce become a staff astronomer at CFHT in 1979 while there he built 
a safer version of the HF cell with sapphire windows for use with the CFHT 
coude spectrograph. The latter was a duplicate of the high dispersion McKel- 
lar spectrograph, and we, together with Stephenson Yang, used the HF sys- 
tem at CFHT over the following twelve years for PRV. The det ector, by then 



an 18 72 diode array, had a correla t ed do uble sampling system (j Walker et al. 



19851 ) based on a design by lGeary! (119791 ) which reduced read-out noise to the 



point where signal to noise was dominated by photon statistics. Full des crip- 
tions of the technique and analyses can be found in lYang et al. and 



Campbell et al.l (119881 ) . To achieve a precision of ~10 m s _1 (equivalent to 
0.006 of a 15/rni diode spacing!) several subtle effects had to be calibrated, 
in particular the impact of convolving the illumination pattern at the slit 
with the HF line profiles. Ultimately, there also turned out to be small but 
significant systematic off-sets in the overall radial velocities between observ- 
ing runs which could be corrected iteratively only at the end of the whole 
program. 

3. The UBC PRV Program 



Yang et al.l (119821 ) carried out the first successful initial program in 1980 
using the HF absorption cell on the CFHT to determine the low amplitude 
velocity variations of the S Scuti star j3 Cass. This was before the high 
reflectance coude mirror train was installed. The wavelengths of the HF lines, 
close to A8700, coincide with a region of reduced reflectivity for aluminum 
resulting in only some 10% overall transmission, a situation soon remedied 
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by the introduction of enhanced silver coatings. This preliminary exercise 
gave us experience and helped to improve the associated reduction programs. 
The results for (3 Cass demonstrated to the satisfaction of successive telescope 
time assignment committees that the HF PRV technique was indeed viable. 

There were two parts to the PRV observing list. Five bright giant and 
one super-giant K stars were included because they were regulary observed 
around the world as radial velocity standards. The basic observing list con- 
sisted of 23 solar m ain sequence and sub-giant stars visible from CFHT 



( jWalker et all Il995l ). Known spectroscopic binaries, rapid rotators and stars 
earlier than F5 were excluded. The stars are listed in Table 1. n is the par- 
allax in asec, the I magnitude corresponds most closely to the wavelength of 
the HF spectral region, and the total number of observations made of each 
star is in the final column. 

We applied twice each year for telescope time and were generally assigned 
four pairs of nights per year - although one year we received none in the first 
six months and the TAC forced us eventually to drop the giants from the 
program and reduced us to three pairs of nights per year. With the mindset 
that we were waiting for planetary signals with characteristic signals of a 
decade, it really was tedious because there could be no obvious results from 
any one observing run and, perhaps more seriously, no publications to nourish 
research funds. Nonetheless, we did persist for 12 years and, it seems, made 
some exo-planetary discoveries in the process - but it took some time for 
these to dawn because none was a simple case. 

3.1. k 1 Ceti and the X8662 Ca II chromo spheric index. 

As mentioned earlier, the chromospheric activity of the individual pro- 
gram stars could be quantified from changes in the equivalent width of the 
chromospherically-sensitive, singly ionised calcium triplet lines one of which, 
at A8662, w as included in our spe ctra. Details of how the index was formed 



are given in lLarson et al.l (jl993bl ). Basically, differences from a mean Ca II 
profile derived from all of the spectra were expressed as AEW 866 . 2 in picome- 
ter units (1 pm = 10 -2 A). 

We expected that there might be a correlation between AEWg66.2 and 
excursions in the PRV curve. Our observations of k 1 Cet appeared to support 
this idea in dramatic fashion as one can see in Figure 2, where AEWs66.2 and 
the PRV values are scaled and plotted on the same diagram. The errors on 
the values of AEWg66.2 are similar to the size of the points while for the PRV 
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Table 1: The UBC PRV Program Stars 
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values they are ~ ±13 m s _1 . The two quantities track each other surprisingly 
well, particularly through the, as yet unexplained, S-shape excursion in 1988. 

In Figure 3 the individual PRV and AEWs66.2 are plotted against each 
other. Despite the increasing scatter at large PRV values, there is clearly a 
significant correlation which was to lead to legitimate doubts about the ap- 
parent planetary PRV signatures for some chromospherically active stars and 
whether a correction based on the k 1 Cet correlation was valid. Ultimately, 
we did not apply any such correction because there was no clear case for any 
other star. 



Later observations of k 1 Cet by iRucinski et al.l (120041 ) and IWalker et al. 



( 120071 ) showed that not only is the star chromospherically active, but that it 
is spotted and rotates differentially and the chromospheric activity is modu- 
lated at a rotational period of 9.3 d. They found no obvious rotational mod- 
ulation of ra dial velocity and point ed out that the <jRy=24.4 m s" 1 found 
over 11 yr by lCumming et al.l (119991 ) would allow a Jovian planet in a tidally 
locked orbit to remain undiscovered. One might also speculate that the sig- 
nificant PRV and AEWg66.2 excursion seen in 1988 lasting some 30 d might 
have been induced by close periastron pas sage of a pla n et in a highly eccen- 
tric orbit by the mec hanism suggested by ICuntz et al.l (120001 ) and for which 
Shkolnik et al.l (120051 ) later found evidence in HD 179949 and v And. 

For good reason, current PRV searches try to avoid chromospherically 
active stars, but that does not mean they do not harbour planets! 

3.2. The Giants 

From initial PRV results one had to conclude that the so-called standards 
were all variable in radial velocity while most of the program stars were con- 
stant! While the range of variability in the giants was not serious enough to 
prevent their use as standards for conventional photographic radial veloci- 
ties, their radial velocities varied between 30 and 300 m s _1 with periods ~1 
year. We published these results ( Walker et al. . fl~98~9l ) after 5 years to point 
out that the yellow giants were a new class of radial velocity variable. In 
particular, Arcturu s showed a range of periods from days to nearly two years 
( llrwin et al.l . Il989l ). Because of the sizes of the giant stars, it was assumed 
that their long rotation periods contributed part of the long-term variability. 



3.2.1. (3 Gem 

As PRV results for the giants accumulated, the repitition of some vari- 
ability bore the signature expected from perturbation by Jovian companions 
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in circular orbits, particularly in the case of f3 Gem (ILarson et all Il993al ) 
where we eventually extended the CFHT PRV series to mid-1993 with HF 
observations from the DAO. For f3 Gem, some eight cycles of a 584.7 day 
period emerged with an amplitude of 46 m s -1 that could correspond to a 
co mpanion with Msim=2.5 M.j (see Figure 4). 



Hatzes &: Cochran! (119931 ) had also reported variable radial velocities for 



three K giants, one of which was Pollux(=/3 Gem) and, in 2006. lHatzes et al 



(120061 ) published the 25 year radial velocity curve shown in Figure 4. The 
velocities include the CFHT and DAO HF PRV values and extended obser- 
vations from the McDonald and Thuringia State, Observatories. There is a 
clear 589.6 day period with an amplitude of 41 m s _1 which, for a 1.7 solar- 
mass primary corresponds to M p ;sim=2.3 Mj with the planet in a nearly 
circular orbit (e = 0.02) and a semi-major axis of 1.6 AU. The persistance 
of the 589.6 day signal ruled out pulsational or rotational variability of the 
star as plausible, alternate explanations. There seems little doubt about the 
reality of this planetary identification. 

3.3. The Main Sequence Stars 



After six years, we ([Campbell et all Il988l ) published interim PRV results 
for 12 solar main sequence stars and 4 sub-giants. Unfortunately, our sug- 
gestion in that paper that statistically significant long-term trends for seven 
of the stars might represent 'the tip of the planetary mass spectrum' only 
seemed to intensify skepticism about the possibility of extra-solar Jovian 
planets. 7 Cep and e Eri were two of the candidates. 

Sometime later, Bruce Campbell, who had been unable to secure a per- 
manent position in astronomy, left the project. Ala n Irwin stepped in to help 



with the data analysis and we (IWalker et all Il995l ) eventually published a fi- 



nal analysis for 22 stars based on the full 12 years of observations. (Published 
in Icarus because they had no page ch arges and by then we had exhausted 
our publication funds!) By that time, Wolszczan fc Frail ( 19921 ) had pub- 
lished their detection of terrestrial-mass planets orbiting the 6.2-ms pulsar 
PSR B1257 + 12, but it was just prior to the dramatic ann ouncement of the 
Jovian -mass planet orbiting 51 Peg with a 4.23 d period by lMayor fc Queloz 
(119951 ). Ironically, the latter discovery was made without the benefit of cap- 
tive gas fiducials but relied on a highly stable spectrograph with a stable line 
spread function and comparison arc lines. 

There were several obvious conclusions from our twelve years of PRV ob- 
servations. While we had detected two previously unknown spectroscopic 
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binaries (7 Cep and x 1 Ori) there was no evidence for brown dwarf com- 
pa nions to any o f our stars similar to, or more massive than, that found 
by lLatham et al.l (119891 ) associated with HD 114762. This apparent brown 
dwar f 'desert' had turned up in other studies (e.g. iDuquennoy fc Mayorl . 
1991c iFischer fc Marcyl . Il992l ). Several stars such as r Cet had shown sur- 
prisingly little change in radial velocity over the 12 years while others had 
proved much more noisy suggesting a higher level of surface activity. 

In our 1995 paper, we struck a more conservative note than in 1988, con- 
cluding that there was no strong evidence for any Jovian planets in circular 
orbits with periods of years associated with any of our stars - a conclusion 
reinforced by the lack of detections by other groups. However, we did draw 
attention to two candidates, 7 Cep and e Eri, where the stellar velocity curves 
could indeed be interpreted as due to Jovian planets. Neither case was sim- 
ple, which is why they remained in doubt until th e excellent follow-up wor k 
of the Texas planet searchers led by Bill Cochran (jWittenmyer et al.l . 120061 ). 



3.4- 7 Cephei 



We had found (jCampbell et al.l . Il988l ) that 7 Cep was a single-line spec- 



troscopic binary from a 'large' systematic change in its radial velocity, but 
our obser vations c overe d only part of the binary orbit (see Figure 5). Sub- 
sequently, iTorred (120071 ) was able to assign an orbital period of 67 yr to the 
system from mu c h mo re extended radial velocity measurements. Recently, 
Neuhauser et al.l (120071 ) actually imaged the faint second star with the aid of 
adaptive optics on the Subaru Telescope - it is likely an M4 dwarf. 

The interest to us was less in the stellar binary than in the 2.5 yr 'ripple' 
supe rimposed on t he ve locity curve seen best perhaps in Figure 5 taken 
from IWalker et al.l (119921 ) that shows the radial velocity curve from the first 
11 years of the program. On removing a low-order-polynomial to account 
for the reflex motion of the primary caused the lon g period stellar binary 
orbit, there was an underlying periodicity from which ICampbell et al.l (119881 ) 
had suggested 'Probable third body variation of 25 m s _1 amplitude, 2.7 yr 
period, superimposed on a large velocity gradient'. With the 11 year curve, 
the periodicity in the residuals was much clearer, and the period was refined 
to 2.52 years. 

But was it a planet or an effect of the star's rotation? At the time, ro- 
tation seemed a slightly more plausible explanation because of a very weak 
almost synchronous chromospheric activity signal coupled with a misclassifi- 
cation of the primary as a giant. There was also the implausibility that there 
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would be a stable planetary orbit in such a close stellar binary, and the then 
prevelant view that Jovian planets would not survive so close to such a large 
primary Unfortunately, o ur best guess at a binary orb i tal pe riod of ~30 
yr later proved too short (ITorresl . 120071 ). and iFuhrmannl (120031 ) established 
that the primary was a Kl IV sub-giant, not a giant. Both of these devel- 
opments undercut any stellar origin for the periodic variations. Ultimately 
Hatzes et al.l (120031 ). by extending measurements of PRV and chromospheric 
activity to 2002 at the McDonald Observatory, demonstrated that the 2.5 
yr period in the residuals had repeated like clockwork and that there was in 
fact no long-term synchronous chromospheric activity. This is shown dramat- 
ically in Figure 6 where the 22 years of differential radial velocities corrected 
for the velocity changes induced by the stellar binary are shown for each 
of the CFHT and three McDonald programs. Zero-point corrections were 
introduced betw een each of th e four series for continuity. 

According to Torres! ( 2007 ). the M4 dwarf secondary star is in an eccentric 
66.8 yr orbit (e = 0.41) with a semimajor axis of 19 AU. The minimum mass 
of the planetary companion to the primary is 1.4 Mj. The planetary orbit 
is only 9.8 times smaller than the orbit of the secondary star (the smallest 
ratio known among exoplanet host stars in multiple systems at this time), 
but it is stable if the two orbits are coplanar. 

So, arguably we detected the first exo-planet PRV signature for a Jo- 
vian planet in a quite unexpected system - a binary star, but it was at a 
characteristically Jovian distance from its primary. 



3.5. e En 



Eri is a young (jSoderblom fc Dappenl . Il989l ). magnetically active star 



(I Gray fc Baliunad . Il995l ) which rotates t wice as fast as the Sun. It also 
has spots and the rotation is differential (jCroll et al.l . 120061 ). As a result, 



its integrated surface velocity changes make it a 'noisy' target producing a 
considerable intrinsic scatter in radial velocity. Despite the 'noise', there 
appeared to be a consistent ~10 yr period in the UBC PRV, but the time 
baseline was too short to confirm this. From observations at McDonald, Lick 



and ESO observatories, lHatzes et al.l (120001 ) extended coverage to 2000. The 
periodicity was still present, and they derived a period of 6.9 yr with an 
amplitude of 19 m s _1 corresponding to a minimum planetary mass of 0.9 
Jupiter in an eccentric orbit (e=0.6). 



Fi gure 7 shows the radial velocity curve extended to 2006 by Benedict et al 
(120061 ). who combined these values with three years of astrometric measure- 
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ments of milliarsecond precision taken with the HST. The noisiness of the 
radial velocity observations is obvious and the astrometry covers less than 
half an orbit. 



4. Conclusions 



By 1993, when our PRV program ended, the HF technique had passed its 
prime. Our observing list was necessarily small because the 100 A coverage 
was an order of magnitude less th an with I 2 vapor and CCD arrays coupled 



to an echelle sp ectral format (e.g . iMarcv fc Butlerl . Il992l ; IWittenmyer et al. 



20061 ). Indeed, iMayor fc Queloa (119951 ) had even demonstrated that with 



sufficient care one could perhaps dispense with an absorption cell altogether. 
The handling and preparation of the HF system was also dangerous. On the 
other hand, the HF system had been straight-forward to implement. The 
natural width and spacing of the HF lines had avoided the additional layer 
of processing necessary with the sharp overlapping I 2 lines. The photometric 
fidelity of the diode arrays was ideal, and we already had a 'small-shift' 
progam in place as a basis for the analysis. Our early success in reducing 
PRV errors had, we liked to think, encouraged other groups to press on. 

As with any new level of precision, our PRV turned up unexpected re- 
sults such as the variability of the yellow giants, the correlated variations of 
the PRV and chromospheric activity in k, 1 Cet, and a paucity of close brown 
dwarf companions. These finds complicated interpretation of our results for 
the program stars. Above all, a general skepticism still greeted any announce- 
ment of an exo-planet discovery. Nonetheless, three stars that we flagged as 
possible planet hosts, 7 Cep, (3 Gem, and e Eri were later confirmed by the 
Texas group, albeit not everyone is yet convinced ab out e Eri b despite th e 
complementary astrometric acceleration detected by iBenedict et al.l (120061 ). 
In simple terms, roughly some 10% of our program stars turned out to have 
Jovian planets with orbital periods of years. This roughly matches the statis- 
tics from the current huge samp le of stars being monitored by others (see for 
example IMarcv fc Butlerl . |2000| ) , although each of our three candidates was 
unusual. 

Since our stars were only observed a few times per year, we had little 
chance of detecting closely bound systems like 51 Peg, but the results for k, 1 
Cet are tantalising in that respect. 

The success of our extra-solar planet search depended on many people, 
but particularly on Bruce Campbell, Stephenson Yang, Alan Irwin and Ron 
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Johnson who built successive low-noise versions of the Reticon camera sys- 
tems. It is a great pleasure to be able to acknowledge my debt to them 
here. 
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Figure 1: Two spectra of 61 Cyg A taken with the CFHT coude spectrograph. The lower 
one was taken with the 80 cm hydrogen fluoride absorption cell in the beam and shows 
the 3-0 HF band lines superimposed, while the upper one was taken without the HF cell. 
The broad stellar feature at A8662 is one of the Ca II triplet lines the strength of w hich 
was used as a measure of chromospheric activity. (Taken from ICampbell et all 119881 ) 
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Figure 2: This double plot demonstrates how the precise radial velocities (PRV) of k 1 Cet 
and its chromospheric activity index (AEW866.2) varied in unison over eleven years. The 
effect is particularly marked through the S-shaped change in 1988. The uncertainties in 
the AEW866.2 values are similar t o the size of the poin ts and the uncertainties in the PRV 



are ~ ±13 m s . (Adapted from lWalker et al.l . ll995f ) 
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Figure 3: The AEW866.2 values for n 1 Cet in Figure 2 plotted against PRV demonstrating 
a significant correlation. 
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Figure 4: Precise radial velocities for the KO giant (3 Gem taken over 25 years from 6 data 
sets, CFHT (green diamonds), DAO (red filled circles), three sets from the McDonald 
Observatory (green open circles, red inverted triangles and bluesquares), and the Thuringia 
State Observatory (red triangles). The solid line corresponds to a 589.6 day period orbit 
and a minimum planetary mass of 2.3 Ju piter and a nearly ci rcular orbit (e = 0.02) with 



a semi-major axis of 1.6 AU. (Taken from Hatzes et all 12006) 
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Figure 5: The 11-year CFHT PRV curve for 7 Cep (|Walker et all Il992h . The overall 



change of ~1 km s 1 is caused by the reflex motion of the p rimary induced by a much 



fainter secondary star (later detected by lNeuhauser et aL . 2007 ) in an orbit of decades. The 



ripple superimposed on the binary star velocity trend was first noted bv lCampbell et al 



(|1988I ) and corresponds to reflex motion of the primary in response to a companion Jupiter- 
mass planet. 
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Figure 6: Residual velocity measurements from four data sets covering 21 years from 
CFHT and McDonald for 7 Cep after subtracting the stellar binary motion. For the 
planet the orbital solution corresponds to a minimum ma s s of 1 .7 Jupiter in a 2.5 yr (2.1 
AU) eccentric orbit (e = 0.12). (Taken from lHatzes et all [20031) 
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Figure 7: The 26 yr radial velocity cur ve for e Eri based on observations from CFHT, 
McDonald, Lick and ESO Observatories ((Benedict et al. . 2006). The star is chromospher- 



ically and magnetically active and consequently a 'noisy' target. A 6.9 yr periodicity 
seems to have persisted with an amplitude of 19 m s _1 which corresponds to a minimum 
planetary mass of 0.9 Jupiter in an eccentric orbit (e=0.6) with a semi-major axis of 3.4 
AU. 
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